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ABSTRACT 

We quantify possible differences between turbulent dynamo action in the Sun and the dynamo 
action studied in idealized simulations. For this purpose we compare Fourier-space shell-to-shell 
energy transfer rates of three incrementally more complex dynamo simulations: an incompress- 
ible, periodic simulation driven by random flow, a simulation of Boussinesq convection, and a 
simulation of fully compressible convection that includes physics relevant to the near-surface 
layers of the Sun. For each of the simulations studied, we find that the dynamo mechanism is 
universal in the kinematic regime because energy is transferred from the turbulent flow to the 
magnetic field from wavenumbers in the inertial range of the energy spectrum. The addition of 
physical effects relevant to the solar near-surface layers, including stratification, compressibility, 
partial ionization, and radiative energy transport, does not appear to affect the nature of the 
dynamo mechanism. The role of inertial-range shear stresses in magnetic field amplification is 
independent from outer-scale circumstances, including forcing and stratification. Although the 
shell-to-shell energy transfer functions have similar properties to those seen in mean-flow driven 
dynamos in each simulation studied, the saturated states of these simulations are not universal 
because the flow at the driving wavenumbers is a significant source of energy for the magnetic 
field. 



Subject headings: magnetic fields 
dynamo 

1. INTRODUCTION 



magnetohydrodynamics (MHD) - Sun: photosphere - turbulence 



A turbulent 3-D flow in an electrically con- 
ducting fluid is capable of generating a magnetic 
field by self-excited dynamo action if the magnetic 
Reynolds number is sufficiently large. In the ab- 
sence of rotation or large-scale shear flow, the en- 
ergy of the generated magnetic field resides pre- 
dominantly at spatial scales significantly smaller 
than the driving (integral) scale of the turbu- 
lence ( [Schekochihin et aL]|2QQ4 ). Such dynamos 
are therefore usually referred to as small-scale dy- 
namos (SSDs). Turbulent plasmas are found in 
nearly all astrophysical systems. Since the big spa- 



tial scales of these systems imply large Reynolds 
numbers, SSDs are expected to operate in most 
of them. This leads to a widespread magnetiza- 
tion of astrophysical plasmas, the consequences of 



which can be profound (e.g.. 


Schekochihin & Cow- 


ley|2006t Wang & Abel|2009t 


Schleicher et al.|2010t 


King & Pringle||2010t Ryu et al.||2008). 



Owing to their small-scale nature, observational 
evidence for magnetic fields generated by SSDs 
is difficult to obtain. Methods have been sug- 
gested to obtain information about the turbulent 
magnetic field in the interstellar and intracluster 
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media (e.g., Waelkens et al. 2009). Evidence for 
such a field with mixed polarity on subresolution 
scales in the solar photosphere is obtained using 
the Hanle e ffect ( |Trujil lo Bueno et al.|2004| |Kleint 
et al.|[2QTQ ) and, possibly, by considering the sta- 
tistical properties of the resolved fields measured 



through the Zeeman effect (Pietarila Graham et al. 
2QQ9a|bD . 
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Apart from these observational attempts, most 
research into SSDs has fohowed theoretical (e.g., 
Kazantsev 1968) or numerical simulation ap- 



proaches. Simulations of SSD action were car- 
ried out for a variety of physical settings, such 
as forced homogeneous incompressible turbulence 



(Meneguzzi et al. 1981j), Boussinesq convection 
( Cattaneo|1999 ), anelastic convection (Brun et al. 
2004), up to fully compressible convection includ- 



ing relevant physics in the solar near-surface lay- 
ers, such as radiative transport and partial ion- 



ization (Vogler & Schiissler 2007 Pietarila Gra- 
ham et al. 2010). Given such a range of physi- 
cal conditions, the question arises of whether the 
mechanism of SSD action is universal or is qual- 
itatively different in the presence of additional 
physics, e.g., as present in the near-surface layers 
of the Sun. In this paper, we therefore analyze 
three incrementally more complex simulations of 
SSDs, namely, (1) a simulation of forced homo- 
geneous incompressible MHD turbulence, (2) a 
simulation of Boussinesq (incompressible) convec- 
tion, and (3) a simulation of compressible and 
stratified solar near-surface convection. To com- 
pare the simulations, we consider energy spectra 
and shell-to-shell energy transfer rates in the kine- 
matic growth phase and in the saturated state of 
the dynamo. Shell-to-shell energy transfer anal- 
yses measure the exchanges of kinetic and mag- 
netic energies between different wavenumbers. In 
the case of incompressible MHD, the method has 
been well studied for dynamos as well as for de- 
cayin g MHD fiows ([Par et al.|2001||Debhquy et al. 



2005 



2010 



Mininni et al.||2005| |Carati et al.||2006| |Cho 



Schekochihin et al. (2007) have raised the ques 



tion whether there is an essential physical differ- 
ence between incompressible mean-flow-driven dy- 
namos dAlexakis et al. 2005[ [Mininni et al.||2005D 
and those driven by random flows with correla- 
tion times shorter than their own turnover times. 
They suggested that, by measuring the shell-to- 
shell transfer of a dynamo resulting from using 
the latter forcing, one should be able to settle this 
question in the following way: if inertial-range 
motions dominate the ampliflcation of the mag- 
netic fleld, the dynamo is purely a property of 
the inertial range and independent of any system- 
dependent outer-scale circumstances. We call this 
condition universal in the Kolmogorov sense. 



Schekochihin et al.|(|2007|) were concerned about 

tions in 
(|2005l); 



the large role of driving-scale motions in the dy- 
namos studied by 



Alexakis et al. 



Mininni 



et al. (2005) and postulated these features were 
peculiar to the mean-flow driven case. How- 
ever, Carati et al. (2006) studied a non- mean-flow- 



driven dynamo for the saturated state and found 
a strong nonlocal contribution from the forcing- 
scale motions to all scales of the magnetic fleld (as 
seen for mean-flow driving). 

The paper is organized as follows. In Sect. [2j 
we describe the three simulation runs and give a 
brief account of the shell-to-shell analysis. This 
method requires a special treatment in the case of 
the solar convection simulation, which is described 
in Appendix A. The results of our analysis are 
reported in Sect. [Sj We discuss them in Sect. [4] 
and give our conclusions in Sect. [5] 

2. METHODS 

2.1. Homogeneous Turbulence (HoT) 



We use a pseudo-spectral FFT code ( Gomez 
et al.||2005a|b| [Mininni et al.||2010[ ) to solve the 
incompressible MHD equations in a periodic box 
with L = 27r, 



dv 

'dt 



- oj X V = —Vp + j X 6 + uW'^v + /, 



da 

'dt 



V • i; = V • a = 0. 



(1) 



Here u? = V x i; is the vorticity. The magnetic 
fleld is given in Alfvenic units, 6 = V x a = 
B / s/pjl^ with the vector potential a calculated in 
the Coulomb gauge. The forcing / is an Ornstein- 
Uhlenbeck process: the amplitudes of the complex 
harmonic modes with 2 < < 3 are evolved in 
time according to 



c -c (a ^i\,f^ 

^k^,ky,k^,t+5t — ^k^,ky,k^,t I J- I +y 

\ '^corr / V '^corr 

(2) 

where Tcorr is the correlation time (taken as unity), 
A is chosen such that Vrms ^1.1, and ^ is normally 
distributed noise with variance 5t. The noise is 
randomized every time step, St = 8 - 10~^. The 
values of the diffusivities are v = r] = 8.8 ■ 10~^ in 
accordance with the resolution of the simulation 
which has 512^ modes (without de-aliasing). The 
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integral scale of the motion is Lq = 1.9, giving 
an integral-scale turnover time of tl = Lo/vrms ^ 
1.7, and (kinetic and magnetic) Reynolds num- 
bers Re = Rem = VrmsLo/v ^ 2400. A non- 
magnetic simulation is run until a turbulent sta- 
tistical steady state is reached after ^Sr^; then, 
a seed field of harmonic modes in the range 20 < 
\k\ < 30 is introduced. After the kinematic phase 
of the magnetic field growth, the run is contin- 
ued (initially with lower resolution and Reynolds 
numbers in order to save computing time) until a 
statistically stationary, saturated state is reached. 

2.2. Boussinesq Convection (BC) 

In the Boussinesq approximation the fiuid is 
treated as incompressible except for the inclusion 
of buoyancy effects related to gravity. The nondi- 
mensionalized equations are 

d(jij 9 

V X (v X + j X 6) = i/V -VO X g, 



dt 



db 

di 



-V X {v xb) = r]V'^h, 



dO 
di 



V-i; = V-6 = 0, Vxh = j. 



(3) 



The temperature represents fiuctuations about 
an equilibrium state with a mean vertical temper- 
ature gradient VTq. The magnetic field is repre- 
sented in Alfvenic units with an Alfven number of 



one. 



The BC simulation is a pseudospectral cal- 
culation performed at a resolution of 512^ in a 
fully periodic box. The amplitude of all modes 
with /c^ = is set to zero to prevent the ex- 
ponential growth of these modes (elevator insta- 
bility, see Calzavarini et al. 2006|) in the verti- 



cally periodic box. This non-restrictive "pseudo- 
Rayleigh-Benard" boundary condition inhibits 
the formation of boundary layers that appear with 
Rayleigh-Benard boundary conditions in a verti- 
cally closed box. 

The BC simulation was carried out for a 
Prandtl number Pr = v / n = 1 and magnetic 
Prandtl number Prm = I'/rj = 2. The magnetic 
Reynolds number, defined in terms of the integral 
scale Lq and Vrnis is Rcm = '^rms^o/^ ^ 4000. 
The Rayleigh number is determined using the 
characteristic length-scale of the vertical temper- 
ature gradient and calculated to he Ra = = 



5.0-10^. Defined in this way, the Rayleigh number 
is not simply comparable to the Rayleigh number 
of a system with defined boundaries, but never- 
theless gives an indication of the balance between 
buoyancy and dissipative forces in the simulation. 

The initial state of the BC simulation consists 
of fully-developed hydrodynamic convection. Ran- 
dom fiuctuations of magnetic field, small com- 
pared to the kinetic energy of the system, are 
seeded into the lowest 16 spectral modes in or- 
der to observe the onset of the linear phase where 
magnetic-field energy grows due to turbulent dy- 
namo action. After nonlinear saturation the sys- 
tem enters an energetically quasi-stationary state. 

2.3. Compressible Solar Convection 

We use results from a dynamo simulation of 
near-surface solar convection carried out with the 
MURaM code dVogler et 'aL]|2005D. The simula- 



tion differs from the more idealized simulations 
of MHD turbulence described above by includ- 
ing physical processes that are relevant for solar 
convection: compressibility, stratification, radia- 
tive energy transport, and partial ionization. The 
equations treated with the MURaM code, written in 
conservation form, are 



dgv 



dB 

~dt 



\B\ 



QVV - 



1 



BB 



2fi I fi 



V • {vB - Bv) = -V X {r]V x B), 



de 
di 



V \ e -\- p 
1 



2/i 



B{v B) 



V -{B xr]V X B)^V -{v-t) 
+V-(i^VT) + ^(^-i;) + Qrad. 



(4) 



where g is the density, p is the gas pressure, and g 
is the gravitational acceleration, vv, BB, vB and 
Bv are dyadic products, and 1 is the 3x3 unit 
matrix. The viscous stress tensor, r, is written 
for a compressible medium with a viscosity coef- 
ficient containing a shock-resolving and a hyper- 
diffusive part (for details, see Vogler et al.||2005" ). 
The total energy density per volume, e, is the sum 
of internal, kinetic and magnetic energy densities. 
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T is the temperature and K the thermal conduc- 
tivity. The source term Qrad, which accounts for 
radiative heating or cooling, is determined by in- 
tegrating the equation of radiative transfer over a 
number of directions for each grid cell as described 
in Vogler et al. (2005). The system of equations 



is completed by the equation of state, which de- 
scribes the relations between the thermodynami- 
cal quantities of a partially ionized fluid. 

The numerical procedure uses centered, fourth- 
order explicit finite differences for the spatial 
derivatives on a uniform Cartesian grid and a 
fourth-order Runge-Kutta scheme for the time 
stepping. The boundary conditions are periodic 
in the horizontal directions with a closed, free- 
slip top boundary and an open lower boundary 
that permits the free in- and outfiow of fiuid (for 
details, see Vogler|[2QQ3l [Vogler et al.||2005 ). 

Here we consider the results of the simulation 
that has previously been presented as Run C in 



Vogler & Schiissler (2007) and Run 2 in Pietar- 



ila Graham et al. (2010). The computational box 



represents a rectangular domain with an extent of 
4.86 Mm in the horizontal directions and 1.4 Mm 
in the vertical direction, covering the range from 
800 km below to 600 km above the optical so- 
lar surfac43 The finite-difference scheme uses 
648 X 648 X 140 grid cells, corresponding to a grid 
spacing of 7.5 x 7.5 x 10 kilometers. 

The magnetic Reynolds number for this simu- 
lation is Re^ = VyuisLo/t] ~ 2000. Because the 
viscous stress tensor r is based on shock-resolving 
and hyperdiffusivity, it is not simple to give ex- 
plicit values of the Reynolds or magnetic Prandtl 
numbers. An estimate of the magnetic Prandtl 
number of between 1 and 2 was derived in I Pietar^ 



ila Graham et al. (2010) by considering the Tay- 
lor scales of the fiow and magnetic field. The 
Reynolds number Re is then in the range of 1000 
to 2000. 

2.4. Shell-to-Shell Transfer Analysis 

This method was first derived for incompress- 



ible Navier-Stokes by Batchelor (1953) and for in- 



compressible MHD by Par et al. (2001). We fol- 



■"^ Unlike the other two simulations presented in this paper, 
the simulation of solar convection is not scale-invariant. 
The equation of state and the radiative transfer depend on 
explicit physical scales. 



low here the exposition by Alexakis et al. (2005) 
who previously applied the analysis to kinematic 
and saturated dynamo states of homogeneous tur- 
bulence driven by a mean fiow. We extend their 
study to random forcing as well as to more phys- 
ically realistic simulations by deriving the com- 
pressible MHD shell-to-shell transfer functions. 

The analysis starts with a decomposition of the 
velocity field and the magnetic field according to 
a{x) = ^^ax(ic), where ax is the part of a 
whose three-dimensional wave vector k in Fourier 
space lies in the range K < \k\ < K-\-l. The inter- 
val (i^, i^+1] is referred to as "shell i^" . Logarith- 
mic binning (7^7^07 7^^^ ^o] for 7 > 1 and inte- 
ger n is required to determine scale-to-scale energy 
transfers ( [Eyink fc Aluie 2009). However, we are 
not seeking to answer questions of scale locality 
of dynamo mechanisms (see, instead, jCarati et al. 
2006). Instead, we seek to determine whether the 



mechanism seen in incompressible homogeneous 
dynamos is also at work in our simulations. 

For this purpose we use linear binning as op- 
posed to the coarser analysis resulting from log- 
arithmic binning; logarithmic binning can be re- 
covered by summing over linear bins. The spectral 
rates of change of the kinetic and magnetic ener- 
gies can then be written as 



d_ 
dt 



dt 



E^{K) = [Tv. + Tbvt + TJv] + ^, (5) 
Q 

Eb {K) = [Tbb + TvBT + Tie] , (6) 



where T = T{Q^ K) are transfer functions: 

rvv(g, ^) = - / ^k{v • v)vQ d^x, (7) 

Tbb{Q.K) = - j BK{v-V)BQd?x, (8) 

rBvT(g, K)= j vk{B • V)Bq d?x, (9) 

rvBT(g, K)= j Bk{B • V)VQ d^x. (10) 

Tiv and Tib represent the effects of dissipative 
heating, where the index "I" stands for internal 
energy. T = T{K) corresponds to the contribu- 
tion of external forces (e.g., gravity). The first 
two indices of the T's denote the energy reservoirs 
involved in the transfer; Trs{Q^K) is the rate of 
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energy transferred from field r in siiell Q to field 
s in shell K. If T^s is positive, energy is received 
by s from r (transfer Q K). If it is negative, 
energy is lost by 5 to r (transfer Q ^ K). The 
third index in TbvT and T^bt denotes the medi- 
ating force, in this case magnetic tension. Since 
Eqs. express energy transfers from field 

r in shell Q to field s in shell they satisfy the 
antisymmetry /conservation relations 



%,{Q,K) = -%,{K,Q). 



(11) 



We now introduce the shell-to-shell transfer 
functions for compressible MHD. In the compress- 
ible case, there is work done by magnetic pres- 
sure, indicated by the index "P". The work done 
by (or against) magnetic pressure cannot be sep- 
arated from energy transferred between different 
wavenumbers inside the magnetic energy reservoir. 
Thus, the above system of shell-to-shell transfer 
functions cannot be used in the compressible case. 
As is consistent for antisymmetric pairings, we 
write instead: 

d 



Ek{K) = ^ [7;v + TbvP + Tbvt + TJv] + ^, 
Q 

(12) 

q^Eb{K) = [TvBP + TvBT + Tm] • (13) 



dt 
d 



In addition to viscous dissipation, Tiv now includes 
the internal energy transferred to the kinetic en- 
ergy reservoir by compression. Pietarila Graham 



et al. (2010) show that this new transfer, T^bp, 



accounts for 5% of the magnetic energy generated 
in the MURaM dynamo. 

The transfer of magnetic energy in shell Q to 
kinetic energy in shell K through the magnetic 
tension force is 



1 

2/1 



{pv)k • ( - • ^ 



d^x, (14) 



where the integral is taken over the analysis do- 
main. The transfer of kinetic energy in shell Q to 
magnetic energy in shell K through the magnetic 



tension force is 



2/i 



{B • V)VQ 



d^x. (15) 



The transfers associated with magnetic pressure 
are 



1 



rBvp(0,i^) = -^ / B 



{VK ■ V)Bc 



Q 



(16) 



B 



K ' 



{VQ ■ V)B 



p 



^B V 



V]B^B{V'Vq) (17) 

{p^)q^ 



d^x. 



The transfers due to different force groupings, here 
magnetic tension (index T) and pressure (index 
P), separately satisfy the conservative antisymme- 
try relation, Eq. (11). The application to non- 



periodic boundaries (MURaM) necessitates a win- 
dowing of the data, described in detail in the Ap- 
pendix. 

3. RESULTS 

3.1. Field Morphology and Energy Spec- 
tra 

In the saturated state, the ratio of total mag- 
netic to total kinetic energy is 0.41 for HoT, 0.36 
for BC, and 0.026 for MURaM. The small-scale dy- 
namo is more efficient in the homogeneous cases 
HoT and BC. In the stratified solar case simu- 
lated with MURaM, the restriction of dynamo action 
to the downflow lanes and the losses due to down- 
flows leaving the computational box limit the level 
of magnetic energy to a few percent of the kinetic 
energy. Figure [l] illustrates the structure of the 
velocity field and the dynamo-generated magnetic 
field in the saturated states of the three simula- 
tions. Shown are maps of the vertical magnetic 
field and velocity components on horizontal cuts. 
In all three cases, the magnetic field exhibits the 
typical folded structures arising from small-scale 
dynamo action, with elongated unipolar features 



5 




Fig. 1. — Maps of the vertical components of magnetic field (upper row) and fluid velocity (lower row) on 
horizontal cuts through each simulation box, taken from snapshots in the saturated dynamo state. The 
two magnetic polarities are shaded in turquoise and brown; white corresponds to zero vertical field. In the 
velocity maps, downflows are shaded red and upflows blue. Note that the arbitrary ratio of driving scale to 
box size is a factor of 3 smaller for MURaM. 



and rapid polarity reversals in the transversal di- 
rection, often on a resistive spatial scale. While 
the field structures are fairly homogeneously dis- 
tributed in the HoT and BC cases, the MURaM sim- 
ulation shows an intermittent structure with ex- 
tended patches almost devoid of field. This struc- 
ture results from the up-down asymmetry of con- 
vection in a stratified medium: the convective up- 
flows expand heavily in the horizontal directions 
and thus expel the magnetic flux. Thus in the 
MURaM simulation, magnetic field generation takes 
place mainly in the vicinity of the narrow turbu- 
lent downflows. The narrowness of these down- 
flows and the size of the convection cells relative 
to the arbitrarily chosen size of the simulation box 
contribute to a smaller appearance of the magnetic 
field structures for MURaM presented in Figure [l] 



The greater separation between the driving 
scale and the box size can be quantified in the 
kinetic and magnetic energy spectra for the simu- 
lations shown in Fig. |2j Vertical lines indicate the 
integral scale for the turbulent motions, 

Lo= J k-^E^{k)dk/j E^{k)dk, (18) 

and the Taylor microscale Av, 

\l = ^= [ E^{k)dkl I k^E,ik)dk. 
{^^) Jo I Jo 

(19) 

These scales signify the approximate beginning 
and end of the inertial range for a hydrodynamic 
or kinematic state. Both scales occur at higher 
wave numbers for MURaM. 
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wavelength [kinematic Lq] wavelength [kinematic Lq] wavelength [kinematic Lq] 

1.0 0.1 1.00 0.10 0.01 1.0 0.1 




Fig. 2. — Kinetic and magnetic energy spectra in the kinematic (upper) and saturated (lower) states of the 
homogeneous turbulence (HoT) simulation (left), the Boussinesq convection (BC) simulation (center), and 
the solar convection (MURaM) simulation (right). 



(Lq/Av)^, which is proportional to the effec- 
tive Reynolds number, is 31.5 in the kinematic 
and 22.5 in the saturated state of the Boussi- 
nesq simulation and 16.5 (10.1) in the HoT simu- 
lation. The corresponding values for MURaM are 
15.7 and 12.7, respectively. The spectra of all 
three dynamos are similar. In Fig. [2| the mag- 
netic energy spectra exhibit power laws with pos- 
itive exponents at small wavenumbers, peaking 
beyond Av in the kinematic state. The omni- 
directional spectra for MURaM are calculated em- 
ploying a Tukey window ( Harris|1978 ) in the verti- 
cal direction. The window corresponds to a high- 
resolution/low-dynamic range apodization. This 
apodization increases spectral leakage from much 
stronger disparate frequencies such as the peak 
of kinetic energy at spatial frequency Lq^. Us- 
ing a high-dynamic-range/low-resolution window 
such as the Nuttall window leads to similar spec- 
tral indices. In the stratified MURaM simulation, the 
anisotropic flow complicates any interpretation of 
the power law exponents. In the saturated state, 
magnetic energy dominates over kinetic energy for 
the largest wavenumbers and the peak of the spec- 
trum moves to lower wavenumbers. 



3.2. Shell-to-Shell Transfer 

In order to compare the dynamo mechanism in 
our simulations, we need to look beyond visual and 
spectral similarities and measure in detail the gen- 
eration of magnetic energy. To do this, we employ 



the shell-to-shell analysis described in Section 2.4 



As an example of how to interpret the trans- 
fer functions T{Q^K)^ Fig. [s] shows the cascade 
of kinetic and magnetic energy towards larger 
wavenumbers in the kinematic phase of the BC 
dynamo. In all of the plots, the transfer func- 
tions have been normalized by the maximum value 
of %BT in the saturated state for each type of 
simulation. The normalization allows straightfor- 
ward comparison of the relative changes in T^bt, 
which is important for understanding the satura- 
tion mechanism of the dynamo. The normaliza- 
tion of all transfer functions for each type of sim- 
ulation is consistent; however the absolute values 
of T^BT depend on the total amount of energy. 
Since the total energies vary between the simu- 
lations, the difference in the absolute magnitudes 
between simulations should not be compared. In- 
stead we focus on qualitative comparisons and rel- 
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Fig. 3. — Shell-to-shell transfer from wavenumber 
Q to wavenumber K in the kinematic state of the 
BC simulation. Top: kinetic to kinetic energy. 
Bottom: magnetic to magnetic energy. 



ative changes, both of which are not affected by 
the normalization. As the values of the different 
transfer functions vary strongly, the 2D plots have 
been constructed to extend over three orders of 
magnitude of positive and negative values. The 
absolute upper limit of the color map corresponds 
to the absolute maximum value of the represented 

r. 

In the top panel of Fig. [Sj the transfer of ki- 
netic energy from one shell to another, 7^v(Q, K), 
is positive (indicated by purple color) for Q < 
i.e., shell K receives energy from slightly smaller 
wavenumbers Q. 7^y{Q^K) is negative (orange) 
foi" Q ^ ^1 shell K loses energy to slightly 



larger wavenumbers Q. Kinetic energy is trans- 
ferred to larger wavenumbers as part of the local 
direct cascade of energy. This is seen also in the 
transfer of magnetic energy from one shell to an- 
other Tbb{Q-,K) (bottom panel). This picture is 
the same for both homogeneous (not shown) and 
Boussinesq convection and for both the kinematic 
and saturated (not shown) states (also compare 
with Fig. 3 of |Mininni et al.||20Q5D . 

The kinematic dynamo mechanism is analyzed 
by measuring Xbt, shown in panels on the left- 
hand side of Fig. [4j This transfer is positive for all 
Q < not just for nearby shells as in the cascade 
process in Fig. [Sj All wavenumbers Q < K of the 
fluid motion drive magnetic energy in a given shell 
K. This is the same as seen for dynamos driven 
by a mean flow (Mininni et al. 2005). Here, we 
see the same kinematic dynamo mechanism in all 
3 simulations. At very large Q > K correspond- 
ing to dissipative wavenumbers, there is less neg- 
ative transfer in the MURaM case than in the other 
two cases, where it is already very weak. We con- 
jecture that this difference is due to the differing 
treatment of viscosity in the codes. 

The saturation of the dynamo results from the 
suppression of the velocity fluctuations at larger 
wavenumber by the back-reaction of the magnetic 
field on the flow through the Lorentz force. In 
the panels on the right-hand side of Figure [4] the 
transfer Q > K is negative (indicated by an or- 



ange color). This generalizes the result of Mininni 
et al.| ( 2005| ) not only to random forcing in ho- 



mogeneous turbulence, but also to incompressible 
convection and to realistic solar-like stratified con- 
vection. The wavenumber range where TvBT < 
coincides with the range at which the magnetic en- 
ergy spectrum dominates the kinetic energy spec- 
trum (cf. Fig. [2|. In MURaM, which has a lower 
saturated ratio £^mag/^kin ^ 2.6% and a larger 
scale separation from the box size, this point oc- 
curs at significantly larger wavenumbers than for 
the other two simulations. 

In the saturated state, there is a transfer of ki- 
netic energy to higher wavenumber {K > Q) mag- 
netic energy similar to that seen in the kinematic 
state. However, the transfer is more localized to- 
wards the small wavenumber Q end of the kinetic 



energy spectrum (Alexakis et al. 2005 Mininni 



et al. 2005 Carati et al. 2006). For example, in 



the HoT simulation, the driving wavenumber for 
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Fig. 6. — Cuts at constant K through the 2D map 
in Figure |4] for the kinematic state of the homo- 
geneous turbulence (HoT) simulation for = 3 
(solid), K = 10 (dashed), = 20 (dot-dashed), 
and K = 30 (dotted). 



the fluid motions, Q ^ 3, is the dominant source 
for transfers to all K in the saturated state. 

Figure |5] presents vertical cuts at constant K 
through the 2D maps of 7^bt(Q,^) shown in 
Figure |4] These slices allow us to quantify the 
changes between the kinematic and the satura- 
tion regimes. For the MURaM dynamo, both the 
shift towards smaller Q and the amount of nega- 
tive transfer for Q ^ is less dramatic than for 
the homogeneous and Boussinesq dynamos. This 
is reasonable because the MURelM dynamo is only 
slightly super-critical and thus the dynamo action 
is an order of magnitude weaker than in the other 
two simulations. If we were able to conduct the 
analysis for much larger i?em, the effect would be- 
come more pronounced. An additional HoT simu- 
lation performed in a regime similar to the MURaM 
dynamo confirms this interpretation. This ad- 
ditional weak HoT simulation was conducted at 
Re = Rem ~ 270 and E^ag = 10% £^kin in the 
saturated state. The MURaM dynamo transfers are 
shown next to similarly shaped transfers from the 
weak HoT in the lower panels of Fig. [5] 

In Fig. [5] a dominant peak at the driving 
wavenumber is absent in the kinematic state. 
Such a peak would be expected for a mean-flow 
driven incompressible dynamo (as seen in Fig. 6 



the transfer from kinetic energy at Q ^ 3 
is no longer dominant for > 3, and the total 
transfer to ~ 3 is dwarfed by the transfers to 
larger K. All three of these results present a view 
of randomly-forced dynamos that differs strongly 
from mean-flow driven dynamos. In particular, 
these results suggest that among small-scale tur- 
bulent dynamos randomly-forced dynamos are 
more local than mean- flow driven dynamos. Tur- 
bulent fluctuations play a greater role, for the 
kinematic phase, relative to integral-scale motions 
in the randomly-forced case than in either mean- 
flow-driven case of Mininni et al. (2005). 



of Mininni et al. 2005). For our randomly- forced 
HoT simulation (Fig. 6]), we find that the total 
transfer to a given wavenumber K increases with 



Work against magnetic pressure becomes in- 
volved in the transport of magnetic energy from 
smaller to larger wave numbers, a phenomenon of- 
ten referred to as the magnetic energy cascade. 
For MURaM this cascade creates an imbalance re- 
sulting in an injection of 5% of the magnetic en- 
ergy generated by the usual dynamo mechanism 
due to the magnetic tension force (Pietarila Gra-] 
ham et al.|2010 ). The shell-to-shell analysis of the 
transfer T^bp is shown in Fig. ^ Net energy is lost 
from small wavenumbers and deposited at larger 
wavenumbers. The transfer is not strongly depen- 
dent on the wavenumber Q of the fluid motions up 
to a break-over point, K ^ 100. Magnetic energy 
at larger scales is expelled by the convection and 
compressed into the downflows. Quantitatively, 
T^BP is stronger for smaller Q, indicating the im- 
portance of motions at the convective scales for 
the flux expulsion process. At larger wavenum- 
bers, magnetic energy is generated by fluid mo- 
tions at smaller wavenumber, e.g., through work 
against the magnetic pressure, and lost to motions 
at larger wavenumbers, such as viscously damped 
magnetosonic waves. 

4. DISCUSSION 

In the three simulations presented here, turbu- 
lent flows are driven by different mechanisms: ran- 
dom forcing (HoT), Boussinesq convection (BC) 
and radiative cooling-driven convection (MURaM). 
The HoT and BC cases are very similar in all as- 
pects. This includes the structure of the flow and 
magnetic field (Fig.jl]), the energy spectra (Fig.[2| 
and the energy transfer spectra (Fig. [4|. The 
similarity is not surprising considering that the 
only essential difference between the simulations 
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is the phase information contained in the large- 
scale driving function. 

The outer-scale appearance of the MURaM simu- 
lation is different from the other two cases (Fig.jl]). 
There is a strong asymmetry between the up- 
flowing and the downflowing plasma, with rela- 
tively smooth upflows and narrow, highly turbu- 
lent downflows. Also, the arbitrary ratio between 
the box size and the driving scale is a factor 3 
larger. In spectral space, this corresponds to a 
shift towards higher wave numbers. The shift af- 
fects the peaks of the kinetic energy (Fig. [2| and 
the kinematic transfer functions (Fig. |5|. Since 



the e-folding time of the magnetic energy in the 
kinematic regime is significantly shorter than the 
turnover time of the granular convection, most of 
the magnetic field appears in the downflow lanes. 

The shell-to-shell energy transfer functions 
(Fig. [4| are similar in the kinematic state, in par- 
ticular transfer from the inertial range wavenum- 
bers is dominant. This suggests that the dynamo 
mechanism, namely the turbulent shear stress of 
the motions in the inertial range, is essentially 
the same in each of the cases studied. Neither 
the short correlation time of random forcing nor 
the additional physics, as present in the Sun, al- 
ter the dynamo mechanism previously discovered 
in incompressible turbulence with a mean flow. 
As a property of the inertial range, independent 
of system-dependent outer-scale circumstances, 
the dynamo mechanism is universal in the Kol- 
mogorov sense. 

In the saturated state, the transfers involve 
much stronger energy transfer from flows at the 
driving wavenumber. This is particularly true 
for the strongly supercritical HoT and BC cases. 
Since the flow at these wavenumbers is determined 
by the large-scale driving, the dynamo action is no 
longer universal in the sense defined above. 

The structure of magnetic fields produced by 
the small-scale dynamo for magnetic Prandtl num- 
ber Prm ^ 1 is different than for Pr^ > 1 
( Schekochihin et al. 2QQ7p . We have shown that 



buoyancy and stratification do not significantly 
alter the small-scale dynamo mechanism. Since 
small-scale dynamo action is known to work for 
^ 1 in incompressible MHD ( Ponty et al. 
2005 Schekochihin et al.|2007 ) and with the same 
shell-to-shell transfer mechanism down to at least 
Prm ~ 0.5 ( Alexakis et al.|20Q7 ), we expect small- 
scale dynamo action to be possible for Pr^ <C 1 
both for MURaM and the Sun. 

5. CONCLUSIONS 

We have compared simulations of solar magne- 
to-convection and local dynamo action with the 
properties of idealized systems in order to evalu- 
ate the robustness and the range of applicability 
of the conclusions drawn from studying the ide- 
alized models. The dynamo has similar shell-to- 
shell energy transfer properties for homogeneous- 
isotropic-incompressible turbulence, Boussinesq 
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convection, and solar conditions that include 
stratification, compressibility, partial ionization 
and radiative energy transport. The results sug- 
gest that the dynamo mechanisms, namely, tur- 
bulent shear stresses acting in the inertial range, 
operate in the same way in each of the cases con- 
sidered. 

For incompressible turbulence we find many 
similarities in the dynamo generated by random 
forcing with a correlation time shorter than its 
turnover time and that resulting from mean-fiow 
driving previously reported ( jAlexakis et al. 2QQ5j 



Mininni et al. 2005). While the signature of 



the dynamo mechanism is the same, the role of 
forcing-wavenumber fiuid motions is diminished 
when random forcing is used for the kinematic 
phase. For the saturated state, injection from 
forcing-wavenumbers remains significant even for 
random forcing (similar to what was found by 
Carati et al.||2006| for yet another type of forcing). 
Basic properties of the turbulent small-scale dy- 
namo process have been thoroughly studied for ho- 
mogeneous, isotropic, triply-periodic simulations; 
these properties carry over to two situations that 
include more complex physics: Boussinesq convec- 
tion and solar surface convection. 

This work has been supported by the Max-Planck 
Society in the framework of the Interinstitutional 
Research Initiative Turbulent transport and ion 
heating, reconnection and electron acceleration in 
solar and fusion plasmas of the MPI for Solar Sys- 
tem Research, Katlenburg-Lindau, and the Insti- 
tute for Plasma Physics, Garching (project MIF- 
IF-A-AERO8047). 
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Fig. 4. — Shell-to-shell transfers for homogeneous turbulence (HoT, top panels), Boussinesq convection (BC, 
middle panels) and solar convection (MURaM, bottom panels). Shown are the transfers of kinetic energy at 
wavenumber shell Q to magnetic energy at shell 7^bt(Q, K)^ in the kinematic (left panels) and saturated 
states (right panels). Within each of the three types of simulations, a given set of arbitrary units is employed. 
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Fig. 5. — Cuts at constant K ^ 1.3/Av through the 2D maps in Figure |4] for the homogeneous turbulence 
(HoT) simulation (top left, i^T = 16), the Boussinesq convection (BC) simulation (top right, = 16), the 
solar convection (MURaM) simulation (bottom left, K = 38) and the weak homogeneous turbulence simulation 
(^mag = 10%£^kin7 ^ 270, bottom right, K = 5). Solid lines correspond to the saturated state (left 
vertical axis) and dashed lines to the kinematic state (right vertical axis). 
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A. Shell-to-Shell Analysis for MURaM 



The use of discrete Fourier transforms is complicated by the non-periodicity of the MURaM data in the 
vertical direction. Ringing effects may seriously taint the transfer functions. We prevent ringing effects in 
this calculation by applying a 50% Tukey window ( Harris |p!978 ) on the data in the vertical direction before 
the transfer analysis. In addition, we zero-pad beyond the extent of the original data to exclude wrap-around 
effects. 

We apply the usual shell filter decomposition a{x) = to velocity momentum pv and the 

magnetic field as described in Sect. 2.4 but with a shell width of 4. When computing the isotropic 
wavenumber we relate all components of k to horizontal wavelengths, i.e., K = 1 corresponds to one 



wave cycle in one of the two axes-aligned horizontal directions. The antisymmetry relations, Eq. (11), are 



satisfied analytically, both for periodic boundary conditions and if a window is applied which tapers off to 
zero at the boundary. In general, a surface integral contribution must be applied. 

To assess the reliability of the transfers computed in this way, we performed a series of tests using the 
results of a low-resolution HoT run. In its original form, the data is fully periodic and the transfer analysis 
is correct by definition. If we "truncate" the data by zeroing half of the box in one direction, the transfer 
analysis results in spurious power at high frequencies in the form of ringing noise, see Fig. [8| If, in addition, 
a window is applied to the data, the results obtained are similar to the original. As is illustrated in Fig. [9j 
the slopes and amplitudes are reasonably well approximated. The tests indicate that the transfer functions 
can be trusted within the scope of the study presented in this paper. 

The disparity between the vertical and horizontal extent of the computational domain corrupts the transfer 
functions at low wavenumbers K^Q <7. These large-scale modes either do not exist in the vertical direction 
or are directly affected by the alteration of the data by the window. K ^ 7 corresponds to the vertical 
extent of the unaltered part of the data. In the plots of the transfers for MURaM (Figs. |4j [sj and[7|), only 
wavenumbers above 8 are shown. 
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Fig. 8. — Transfer function in a 128^ HoT simulation computed in three different ways: using the original 
periodic cube (left), using a cube in which half of the data has been zeroed (middle) and using a "half-zeroed" 
cube with the remaining data being forced to be periodic by a Tukey window (right). The results have been 
renormalized to account for the impact of the truncation and windowing on the amplitude. 
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